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Abstract 

We theoretically propose directional, low-threshold plasmon lasing in both the near-infrared and 
visible wavelengths by utilizing spatially coherent surface plasmon polaritons on a meta-surface. 
The gain strength required for threshold lasing can be tuned down to a large extent through 
compatible structural parameters. Our calculations show that no more than 65 cm -1 at 193.5 THz 
(1.55 urn) or 267 cm _1 at 474THz (0.633 /im) of gain coefficient is sufficient to compensate for the 
dissipation of metal films for threshold lasing; these values are smaller than any reported studies 
at the same frequencies. These findings present a planar solid-state route for plasmon lasing that 
is highly efficient and spatially coherent. 
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I. INTRODUCTION 



In recent years, nanophotonics has attracted considerable attention as one of the emerging 
fields of modern science and technology. Active devices, especially coherent sources with high 
efficiency, are of great importance for the development of nanophotonics. The amplification 
of the evanescent field and strong confinement of light in a subwavelength scale are crucial 
for active nanophotonic devices to achieve superior performance. Surface plasmon polaritons 
(SPPs) provide a feasible route for manipulating light beyond the diffraction limit. In 2003, 
Bergman and Stockman proposed surface plasmon amplification by stimulated emission of 
radiation (SPASER) 1,2 . A spaser — referring to the accumulation of a large number of SPPs 
in a single mode assisted by a gain medium — is precisely the plasmonic counterpart of the 
nano-laser. In principle, one of the most important applications of SPP amplification is 
plasmon lasing- 1 ^. As energy dissipation in metals significantly reduces the lifetime of SPPs 
in the infrared and visible frequencies, the gain of the active medium must prevail over the 
dissipation of metals in order to sustain the survival of a spaser. Many efforts related to 
SPP amplification have been devoted to loss compensation in metal-dielectric structures^-. 
The first experimental verification of the spaser used gold nano-particles encapsulated in 
dye-doped silica shells, which support a collective SPP mode with a high quality factor 
(Q-factor)2. 

Planar design is preferred for the realization of spasing as well as plasmon lasing for the 
sake of fabrication convenience and device integration. However, as a finite-sized planar 
structure is usually incapable of strong light confinement, a high-gain medium is required in 
order to compensate for the radiation and damping losses from the planar structure. One 
solution to this difficulty is to utilize "trapped mode" of a plasmonic meta-surface 10 . A split 
ring resonator with weak asymmetry (ASR) possesses a leaky mode that is poorly coupled 
with free-space photons. A plasmonic metamaterial slab composed of an array of such ASRs 
will facilitate SPP amplification and lasing with low-gain material when the leaky mode is 
excited. Subsequent experimental investigation showed that the loss from the nano-sized 
metallic ASRs could be compensated for by the gain of quantum dots^i. 

It is worth noting that a magnetic meta-surface can support spatially coherent SPPs which 
are of high Q-factor as well. The subtle concept, in contrast to the "trapped mode" which 
is derived from dark mode of a single resonator, comes from the SPPs in harmonic modei 2 -. 
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These spatially coherent SPPs are actually the leaky waves of a slab waveguide which can 
not be determined by a single resonator. They must be in harmonic mode to sustain a leaky 
wave character, giving rise to a high Q-factor which can be of one or two orders larger as 
compared to that of the SPPs in fundamental mode. In this paper, we theoretically propose 
that such spatially coherent SPPs present a practical solution for low-threshold plasmon 
lasing in both the near-infrared and visible regimes. The gain coefficient a of the active 
medium for threshold lasing can be tuned down by structural parameters to a large extent, 
even down to values much smaller than those of semiconductor quantum wells. Two model 
samples are presented with structural parameters that are feasible for practical fabrication. 
The threshold gain coefficient a t h =65 cm -1 at 193.5 THz (1.55/zm) is much smaller than 
those of other lasing systems reported in the same infrared frequency region^^^. The low 
threshold arises from the high Q-factor of the SPPs in harmonic mode and from the strong 
confinement of the local field in the active layer, both of which are not easy to establish in a 
planar solid-state system. The high directionality of the lasing beam with the full width at 
half maximum (FWHM) of AO = 4° in E-plane (<p = 0°) and A6 = 16° in H-plane (<p = 90°) 
is verified by finite-difference-in-time-domain (FDTD) simulations, where 9 and ip are the 
polar angle and azimuthal angle. 



II. REFLECTANCE AMPLIFICATION ASSISTED BY SPATIALLY COHERENT 
SPPS 

Our planar lasing system is schematically illustrated in the inset of Fig. 1(a). Lying 
on the xy plane, the slab comprises an upper layer of a metallic lamellar grating with a 
thickness of t, a metallic ground plane, and an active layer with a thickness of h sandwiched 
between them. The metallic strips with a width of a are separated by a small air gap g, 
giving rise to a period of p = a + g for the lamellar grating. The geometric parameters of 
our model in the infrared region are t = 0.05yum, h = 0.2/rni, a = 1.1/rni, g = 0.05/xm, and 
p = a + g = 1.15/im. 

The SPP resonances on the slab can be characterized by the amplification of reflection 
coefficients. The metals in our system are assumed to be silver with permittivity defined 
by a Drude modet^, e m {u) = — Up/(u 2 + icuj) where = 3.7, u p = 13673THz, 
7 = 27.35THz. And the dielectric layer is assumed to be active with a complex permittivity 
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i = e' + is", in which the imaginary part e" is positive. We perform a rigorous coupled- 
wave analysis^ 1 ^ to calculate the coefficient r m (9, /) of the m th order of the reflected Bloch 
wave components, in which 9 and / denote the incident angle and the frequency of plane 
wave incidence. The incident plane wave is transverse magnetic (TM) polarized with the 
magnetic field parallel to the y direction and in-plane wave vector k\\ = kosin9e x , where 
ko is the wave vector in a vacuum. The gain strength of the active layer is defined by a 
coefficient, a = {2n / \)Im(y/ e' + is"). In the calculations, different values of gain strength 
are taken into account by varying the imaginary part e" of the complex permittivity while 
keeping the real part fixed at e' — 1.7. 

Figure 1 shows the calculated reflection spectra in the th order under an incident angle of 
9 = 0°, 10° and an imaginary part of the permittivity of the active layer of e" = 0.878, 0.021. 
When the gain strength is strong {s" = 0.878), there exists an amplification peak fixed at 
46.3 THz (A = 6.48/xm) that is almost independent of the incident angle. In contrast, when 
the gain strength is weak (e" = 0.021), the amplification peak at 46.3 THz disappears, and 
a new peak emerges at 193.5 THz (A = 1.55/im) under normal incidence. Under oblique 
incidence, the emerging peak splits into two peaks with varying frequencies that are sensitive 
to the incident angle; see, for example, the dotted line in Fig. 1(b) for the two peaks at 166.1 
THz and 211.9 THz with 9 = 10°. The angle-dependent amplification peaks are obviously 

10' 



>>10'l 

.5 io'i 
■*-> 

o i 

at .„7 



10 s ' 

50 100 150 200 

Frequency(THz) 





FIG. 1. Reflection spectra under TM polarized plane wave incidence at = 0° (red solid lines) and 10° 
(blue dashed lines) with different values for permittivity of the active medium: (a) e = 1.7 + 0.878i, 
(b) e = 1.7 + 0.021i. The inset in Fig. 1(a) shows the schematic geometry of our sample model. 
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FIG. 2. Reflection amplification at 193.5 THz and 46.3 THz as a function of gain strength of the 
active layer. 




different from the angle-independent peak in regards to the physics of their origin. Their 
time-reversed counterparts are the angle- dependent and angle-independent absorption peaks 
that we have investigated before^. An angle- independent peak, requiring a sufficiently large 
gain/absorption coefficient of the dielectric layer, comes from the excitation of a localized 
SPP resonance in a fundamental mode that is dominated by the th order guided Bloch 
mode in the dielectric. The angle- dependent feature comes from the spatially coherent SPP 
resonance in a harmonic mode that is dictated primarily by the ±l st orders of the guided 
Bloch mode in the dielectric. 

The reflection amplification of the two resonant peaks as a function of the gain strength 
in the active layer is plotted in Fig. 2. It should be noted that a rather low gain strength 
ath = 65cm _1 is sufficient to compensate for the loss of metals and to sustain the amplification 
of the guided waves in the active layer as well as the reflectance amplification at 193.5 THz. 
The amplified reflectance increases exponentially with the increase of a, until a maximal 
amplification of about 1.4 x 10 5 is reached at an optimal value of a opt = 342cm -1 . In 
contrast, we see from Fig. 2 that the reflectance amplification at 46.3 THz requires much 
stronger gain strength, which is unattainable for most of the existing active media at this 
frequency. 

Additional calculations indicate that the gain strength required for the reflection ampli- 
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FIG. 3. (a) The time dependence of the magnetic fields of reflection. Spectra are acquired through 
Fourier transformation of the time-domain results with different time spans (b) t = 0.083 ~ 2ps 
and (c) t = 0.167 ~ 2ps. 

fication at 193.5 THz can be reduced further, provided that a smaller air gap or a thinner 
dielectric layer which matches the gain strength is adopted. The property provides us with 
a good opportunity for finding a balance between ease of fabrication and attainable gain 
strength of the active layer. Usually thermal conductivity of active medium becomes larger 
with respect to higher temperature. A weaker gain strength required for threshold lasing 
also implies that we can find an appropriate structure for plasmon lasing that is robust 
against thermal noise, or intuitively that the thermal noise is somewhat "suppressed" by 
spatially coherent SPPs. This is advantageous for room-temperature plasmon lasing as one 
goal being pursued. At this stage, it is imperative for us to apply the spatially coherent SPP 
modes for low-threshold plasmon lasing as such a SPP mode with a long lifetime, which 
is compatible to the structural parameters, will create extensive adaptability in material 
selection and feasible way for practical implementation of plasmon lasing. 
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FIG. 4. (a) Energy level diagram for the Ytterbium-erbium co-doped material Er:Yb:YCOB. (b) 
Time dependence of the averaged H y over all the spatial points of a xz section of the cavity, (c) 
the magnetic field H y distribution at a typical time (t =0.96ps). 

III. RATE EQUATION ANALYSIS: LASING THRESHOLD AND DYNAMIC 
PROCESSING 

Now we investigate the threshold condition of our lasing system by using the realistic 
material parameters. A widely used gain medium, Ytterbium-erbium co-doped material 
Er:Yb:YCOB is adopted as the active medium for lasing around 1.5/im^ 1 ^. Its level assign- 
ment is shown in Fig4 (a). The population density on every level evolves according to the 
dynamic rate equations 20 , and can be solved for steady-state conditions: 

dN 2Y _ dN 1Y 

dt dt U 

Noy 

= R P k FT N 2Y N 1E + k BT N 1Y N 3E = 0, 
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(3) 



dN 3E 



= k FT N 2Y N 1E - k BT N 1Y N 3E + /3 43 



(4) 



dN AE 
dt 




N AE 



(5) 



N ix and t, are the population density and radiative lifetime of corresponding levels, are 
the branching ratio for transitions from the i th to j th levels. R p is the pump rate from 1Y to 
2Y level, &ft and /cbt are the forward and backward energy transfer coefficient corresponding 
to 2Y+1E— >3E+1Y and backwards. C up is the coefficient of the cooperative upconversion 
2E+2E-*4E+1E. W 2 \ is the stimulated-emission coefficient. AN 2 i = iV 2 E - N XE is the 
inverted population difference between IE and 2E level. Ni is the population density of the 
lasing mode. Considering the conservation of population density, we have Nyb — N 1Y + -^2Y 
and N Er = N 1E + N 2E + N 3E + N iE . The threshold lasing can be evaluated under steady-state 
condition as ^ = (W 21 AN th — ^)N t = 0, where Q is the Q-factor of the lasing mode. We 
note that the threshold inversion density AiV t h can be estimated by Q; and a higher Q-factor 
as a result of lower radiation loss and damping loss of the lasing mode will give rise to lower 
threshold doping concentrations or pumping power dominated by AiV t h- 

The Q-factor of SPP mode can be extracted from the FDTD numerical simulations by 
assuming that the dielectric is passive and lossless with a dielectric constant equal to 1.7, 
the real part of that of Er:Yb:YCOB. A one-way TM polarized plane wave with a Gaussian 
distribution in time domain is normally incident on the model slab. The incident plane is 
2.5/xm away from the top surface of the model slab. A probing plane is positioned 0.35/im 
behind the incident plane to monitor the reflected wave fields. Figure 3(a) illustrates the 
time-domain spectroscopy of the magnetic local field H y at the probing plane. The frequency 
spectra in Figs. 3(b) and 3(c) present the Fourier transformed data of the time-domain 
spectroscopy in different time spans t = 0.083 ~ 2ps and t = 0.167 ~ 2ps. A broad peak 
and a narrow one, centered at 45.0THz and 192.6THz respectively, are observed in both 
figures. More calculations in time domain indicate that the decay rate of the broad one is 
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much quicker than that of the narrow one. These two peaks correspond to the localized 
SPP resonance in fundamental mode at 46.3THz and the spatially coherent SPP resonance 
in harmonic mode at 193.5THz, as shown in Figs. 1 and 2. Q-factor of a SPP resonance 
can be extracted from the FWHM in the frequency spectra. A low Q-factor of only 4.2 and 
a high Q-factor of 94.2 are derived at 45.0THz and 192.6THz in the frequency spectra. The 
Q-factor can also be estimated with the exponential decay rate of the mode in time domain, 
as Q = 2irfT, where r is the lifetime of the mode. In our case, the life time of each mode is, 
r = 0.0149ps and 0.0779ps, respectively, reproducing exactly the same values of Q-factor at 
the two resonance frequencies as the former method. We note that the Q-factor of spatially 
coherent SPP mode is of one order larger than that of localized SPP mode. And a higher 
Q-factor will be advantageous in the field amplification in gain medium and the stimulated 
emission radiations. By applying the condition of threshold lasing at 192.6THz to the rate 
equations, a threshold doping concentration 1.48 x 10 25 ions/m 3 of Er 3+ is obtained, which 
falls inside the typical concentration range of Er 3+ ions in Er:Yb:YCOB. In our calculations, 
the Yb 3+ concentration is assumed to be fixed at 5.0 x 10 27 ions/m 3 . It is not likely to utilize 
a fundamental SPP mode of the meta-surface to realize plasmon lasing due to the low Q- 
factor of 4.2, as the required Er 3+ concentration 1.3 x 10 27 ions/m 3 is hard to achieve in 
Er:Yb:YCOB. More calculations show that it is not very easy for the Q-factor of the local 
SPP resonance in fundamental mode to exceed a value of 10 by tuning structural parameters. 

Now we perform FDTD simulations to dynamically demonstrate a lasing example of our 
model slab 21 =~. Periodic boundary is adopted in our simulations, the Yb 3+ concentration 
and the Er 3+ concentration are chosen at 5.0xl0 27 irons/m 3 and 4.0xl0 25 irons/m 3 , respec- 
tively. Calculations are performed by assuming a pumping power at P pump = 2.1 x 10~ 5 W 
which is beyond the threshold P t h = 1-9 x 10~ 5 W . The magnetic field H y pattern of the 
cavity and its emission field distribution at a typical time (t = 0.96ps) are shown in Fig. 
4(c). Fig. 4(b) shows the time dependence of the averaged H y over all the spatial points in 
xz plane covering the unit cell. 

The lasing at 192.6THz of a real structure is also considered. In the FDTD simulations, 
the lateral size of the model slab is 20 periods along x direction (the direction of the grating) 
w x = 20p, and the length of the stripe along y direction is w y = w x . The angular dependence 
of far-field patterns are calculated at the steady-state. We see from Fig. 5 that most of the 
energy of the lasing beam is confined in a small solid angle. The angle-dependent power is 
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FIG. 5. Far-field radiation patterns \Ef at 192.6THz as a function of polar angle 9 and azimuthal 
angle (p. (a) \E(6,tp)\ 2 in a stereogram chart, (b) \E(9,0°)\ 2 (E-plane), (c) \E(6, 90°)| 2 (H-plane). 

normalized to the maximum at the normal direction. The FWHM in the E-plane (<p = 0°) 
is AO = 4°, and that in the H-plane (ip = 90°) is AO = 16°. The highly directional lasing 
beam shall be attributed to the excitation of a spatially coherent SPP mode on the slab, as 
such an SPP mode will inherently establish phase correlations among the unit cells on the 
meta-surface due to the dominant ±l st orders of guided Bloch modes in the active layer. 



IV. LASING IN THE VISIBLE REGION 



Our lasing system can also operate at visible wavelength with an appropriate design. Fig- 
ure 6 presents an example of reflectance amplification at 474 THz (A = 633nm). A rather low 
gain coefficient of a t h = 267cm -1 is enough for the compensation of metal dissipation. Even 
the optimal gain coefficient, a opt = 857cm -1 for the maximum reflectance, falls in the range 
of the gain coefficients of the most commonly used active media in this frequency regime, 
such as organic dye molecules, semiconductor quantum wells, and nano-crystal quantum 
dots. The structural parameters are as follows: grating period p = 400nm, metallic strip 
width a = 385nm, air gap width g = 15nm, grating thickness t = 20nm, and dielectric 
layer thickness h = lOOnm. This calculation demonstrates that the spatial coherence of a 
harmonic magnetic SPP resonance can be harvested for plasmon lasing in the visible region. 
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FIG. 6. The spectra of reflection amplification around 474THz with respect to different values of 
gain coefficient a of the active medium. 



V. CONCLUSION 



In summary, we demonstrate that spatially coherent SPP modes in harmonic mode on a 
planar plasmonic meta-surface can be utilized for low-threshold directional plasmon lasing. 
As both the Q-factor of the spatially coherent SPP modes and the confinement of the local 
field are dictated primarily by structural parameters, the SPP amplification can be triggered 
at a rather low threshold gain strength attainable with most of the active media. Meanwhile 
our examples indicate that the planar low-threshold plasmon lasing is not difficult to be 
realized with current fabrication techniques, even in the visible region. Our calculations 
show that using the harmonic SPP modes is crux of the matter for achieving low-threshold 
lasing on a magnetic meta-surface rather than using the fundamental ones. 
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